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Abstract 

We compute the ma dependence of lattice renormalization factors for sea 
quark matrix elements. The results differ from the (1 + ma) correction factor 
commonly used for valence quarks and connected current insertions. We find 
that for sea quarks, the correction factors are in general larger and depend 
strongly on the Lorentz structure of the current under consideration. Results 
are presented both for the Wilson action and for the 2-link improved action of 
Hamber and Wu. Phenomenological implications are also briefly discussed in 
two examples. 



1 Introduction 



Many interesting problems in hadronic physics (e.g. the proton spin crisis, the U(l) 
problem, the strange quark content of the nucleon, ...) require the computation of 
a quark current inside of a disconnected^ fermionic loop. Matrix elements in the 
nucleon for example generally involve two types of contribution associated with the 
connected (fig.l.a) or disconnected (fig.l.b) insertion of the quark current. Simi- 
larly, the propagation of flavor-singlet mesons involves, besides the usual connected 
contribution (fig.l.c), a disconnected contribution (fig.l.d). Up to recently, only the 
contributions associated with fig.l.a and c (which we will call "valence" for shortjj) 
were considered in lattice simulations. Contributions from fig.l.b and d (which we 
will call "sea" for short), on the other hand, were considered too difficult to compute. 
Recent developments in the technology of lattice gauge theory have however made 
their computation feasible. Among the techniques currently being used are the Z 2 
noise method [Q, [J and the non-gauge- fixed volume source technique ||, |5| . Progress 
is fast in this area and numerical simulations based on these methods are beginning 
to produce statistically significant results. It therefore becomes essential to compute 
the renormalization factor associated with sea quark matrix elements in order to turn 
the raw numbers from Monte-Carlo simulations into physical quantities. The lattice 
renormalization of bilinear quark operators (in the limit ma — > 0) was considered a 
long time ago ||, but the finite ma corrections were only introduced more recently @ 
along with the tadpole improvement [BJ. Here we would like to get a closer look at the 
finite ma corrections in the particular case of sea quark matrix elements. We expect 
this to be important for the problem at hand. One would like for example to be able 
to measure various strange quark matrix elements in the nucleon. In current simula- 
tions, ma is about 0.1 for a strange quark and therefore finite ma corrections might be 
significant in this case, and indeed they are as we will see below. At the heart of the 
present study is the realization that the finite ma correction factor commonly used 
for valence quarks and connected current insertions is not adequate for the case of sea 
quarks or disconnected current insertions. In section 2, we briefly review the usual 
approach and discuss its limitations. We also make some general comments which 
will be useful in understanding the results presented later. In section 3, we motivate 
and introduce the method that we use to compute the finite ma correction factors for 
sea quarks. The method in itself is rather straightforward and mostly consists in com- 
puting the perturbative response of a quark current to an external gauge field both 
on the lattice and in the continuum. By comparing the two results, the appropriate 
correction factor is then extracted. We have carried out the procedure for the scalar, 
pseudoscalar and axial currents which are of most direct phenomenological relevance. 
The computation of the relevant Feynman diagrams in the continuum is presented 
in section 4, whereas the lattice results are to be found in section 5. The finite ma 

*Here the terms connected or disconnected refer only to the topological classification of the quark 
lines (regardless of the presence of gluonic lines). 

' Actually, the connected insertions contain both valence quarks and cloud quarks and antiquarks 

§• 
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correction factors are then computed and the results are summarized in table 1 for 
the Wilson action and in table 2 for the 2-link improved action of Hamber and Wu 
0. The application of these correction factors to Monte-Carlo data is discussed in 
section 6 for two examples (the computation of the mass of the rf and the determi- 
nation of the strange quark contribution to the mass of the nucleon). In section 7, 
we add some final remarks and briefly discuss further developments. Our conventions 
and Feynman rules are presented in the appendix. 

2 Wilson fermions at finite ma 

It has been customary to estimate the finite ma corrections by comparing the lattice 
and continuum propagators of a free quark at p = [{7| . The result of this computation 
gives a (1 + ma) finite renormalization factor for quark bilinears, where ma is the 
mass appearing in the Lagrangian or more generally the mass obtained after tadpole 
improvement: 

1 K 

ma = — — 4 where k = kUq ~ (1) 

2k 8k c 

The assumptions going into this calculation appear to be justified for valence quarks. 
If the quark is heavy, its 3-momentum is small and it is "almost on mass-shell", 
whereas when ma goes to the correction factor goes back to 1. So the approximation 
is reasonable for a quark propagating forward in time inside a meson or baryon at rest 
(although there seems to be room for improvement at small and intermediate values 
of ma). For sea quark matrix elements, however, the situation is completely different. 
Configuration by configuration, one is computing a closed quark loop in a background 
gauge field (see fig.l.b and d). This involves an integration over the 4-momentum 
of the quark running around the loop. Therefore, the kinematical restrictions (i.e. 
low 3-momentum and proximity to the mass-shell) which were justified for valence 
quarks no longer apply for sea quarks. A more general method, avoiding particular 
kinematical assumptions, will thus be necessary to compute the finite ma corrections 
in this case. 

Before introducing our technique for dealing with sea quarks in section 3, we would 
like to make here some qualitative remarks on the nature of the problem. As is well 
known, the r = theory (naive action) contains 16 fermions instead of 1 and the 
role of the Wilson term is to lift the degeneracy between the various species. In the 
limit ma — > 0, it does so no matter what the value of 0) is. However when ma is 
non-zero, the effectiveness of the Wilson term in lifting the degeneracy becomes more 
limited and r dependent. In this sense, at finite ma, one is dealing with an effective 
theory with a more complicated flavor content. There are 16 species with masses [ |T0|| : 

ma ma + 2r ma + 4r ma + 6r ma + 8r 

{1} {4} {6} {4} {1} 

where the degeneracy of each level is indicated between braces. It is clear that when 
ma becomes of the same order as r, the doublers are no longer much heavier than 
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the "fundamental" fermion. This causes strong lattice artifacts, in the computation 
of closed quark loops and requires the use of correction factors in order to map the 
effective "multiflavor" theory onto QCD. It is precisely the computation of these finite 
renormalizations that we will consider in this paper. We will see that they depend 
on the Lorentz structure of the current under consideration. This is to be expected. 
Take naive fermions for example, matrix elements of the scalar current in a smooth 
external gluonic field would be 16 times what they should be in QCD, whereas matrix 
elements of the axial current would be in the chiral limit because the chiralities of 



the 16 doublers add to |10fl . As we will see below similar trends subsist for Wilson 
fermions at finite ma: matrix elements of the scalar current are overestimated whereas 
those of the axial current are underestimated. That is, at finite ma, the Wilson term 
only takes care of the doubling problem partially and this affects the measurements 
of sea quark matrix elements. 



3 Renormalization of sea quark matrix elements 

Since configuration by configuration, the computation of a sea quark matrix element 
reduces to the computation of an "induced quark current" in a given background 
gauge field, we are led to study the "response functions" to an external gluonic field. 
By computing the matrix element of some current in a given background gauge field 
on the lattice (at finite ma) and comparing with the induced current in the "same" 
external field in the continuum, we can estimate the lattice artifacts at this particular 
value of ma and correct for them. A particularly interesting example was considered 
by Smit and Vink several years ago ||11|| . They investigated the case of the flavor- 
singlet pseudoscalar current where a topological relation gives an exact expression for 
the response function. Indeed, by making use of the anomalous Ward identity 

df+j^ - 2mj 5 = 2n f q (2) 

where q is the topological charge density defined as 

Q = ^^ e nvpaTr{G fu/G pa ) (3) 

and integrating over 4-volume, one obtains a direct relation between the zero momen- 
tum matrix element of j§ (Q5 = J d 4 xj^) and the topological charge (Q t = J d A xq) of 
the gauge field configuration: 

m Q 5 = -n f Q t (4) 

It is then sufficient to take a smooth background field configuration of given topo- 
logical charge on the lattice, to measure the matrix element of Q5 at various values 
of ma and to compare with (4) to be able to deduce the appropriate ma dependent 
correction factor. We will call this factor Kp(ma) [P for pseudoscalar]. Smit and 
Vink then went on to compute the topological susceptibility through correlations of 
< J5(ic)j5(0) >. Taking into account the correction factor, one obtains in terms of 
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the lattice currents: 



X =< Qt > /V = Kp(ma) J (fx < j 5 (x)j 5 (0) > m 2 /nj (5) 

Strictly speaking, one has to be careful at this point, since the correction factor could 
be gauge configuration dependent, in which case normalizing on a single configuration, 
as was done above, would not be adequate. However, for the present computation, 
a consistency check can be used to show that finite ma artifacts are to a good ap- 
proximation configuration independent. We know that the topological susceptibility 
(being a pure glue quantity in the quenched approximation) must be independent of 
the mass of the quarks. An evaluation of (5) on a statistical sample of gauge config- 
urations would however not satisfy this constraint if Kp(ma) wasn't the appropriate 
correction factor for every "important" configuration. The fact that Smit and Vink 
indeed obtain for (5) a mass independent result in their Monte-Carlo simulation there- 
fore can be used as an indication that a single, configuration independent, correction 
factor [Kp(ma)] is sufficient. 

Now we would like to generalize the computation of finite ma corrections to other 
currents. The pseudoscalar current was special because of the existence of an exact 
"topological" relation (4); other cases won't be that straightforward. However, if we 
assume that the ma dependence of the correction factor is the same for all configu- 
rations (as seems to be the case from the above discussion), any background gauge 
field can be used for its evaluation. In particular, we can consider very weak exter- 
nal fields, so that perturbation theory becomes a reliable tool, and compute in this 
manner the relation between lattice and continuum theory. This is the strategy that 
we will adopt here. It turns out that the lowest order perturbative diagrams for the 
expectation value of a quark current ipTip in an external field will be triangle diagrams 
like the one represented in fig.2.a . In this paper, we will restrict our attention to 
matrix elements of zero momentum currents (so that q = —p in fig.2). This simplifies 
the computations and is sufficient for the applications currently envisaged ( "forward" 
matrix elements in the nucleon, for example). 

Let's consider the scalar current as an example (r = 1 in fig.2) . It is well known 
that in this case, the continuum computation of the triangle diagram, gives a relation 
between the matrix element of the scalar current and the square of the external field 
strength of momentum p (the zero momentum relation has been abundantly used in 
the context of QCD sum rules |L2]): 



P 



2 



m<W>>= C(£a)^(p) (6) 



m 



where C is a dimensionless function of p 2 /m 2 . Similarly, one can compute the triangle 
diagram on the lattice (see section 5 for details) and obtain the corresponding relation: 

m < ijjif) >= L(pa,ma)F 2 u (p) (7) 

the only difference being that the lattice dimensionless function L can depend sepa- 
rately on pa and ma. The functions L and C contain all the information we need, 
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they tell us how the lattice and continuum scalar currents react to a given external 
field. By comparing the two, we can deduce which correction factor need to be in- 
troduced in lattice measurements of the scalar current. We will mostly consider the 
case p=0 (corresponding to smooth background gauge fields). Then one obtains the 
finite (ma dependent) renormalization factor for the scalar current: 



Of course, one can check that L(pa, ma) goes back to C(p 2 /m 2 ) when both pa and ma 
are small so that ^(ma) — ► 1 when ma —>■ 0. If the numerical simulation is dominated 
by smooth gauge fields, then Kg (ma) is all that we need. However, this is not in 
general the case and it is therefore important to look also at the situation at finite pa. 
To investigate this momentum dependence, we will present below comparative plots 
of L(pa,ma) and C ((pa) 2 / (ma) 2 ) at fixed values of ma. The perturbative method 
presented here for the scalar current generalizes, mutatis mutandis, to other currents. 
One can for example compute the finite ma correction factor for the pseudoscalar 
current and then compare with the topological charge method presented earlier. The 
two results agree provided that the "instanton" field is weak enough. 

A general remark is in order here. In the case of valence quarks [0, it is assumed 
that the ma and /3(or a s ) dependence of the renormalization factor factorize: 



We make the same factorization hypothesis here. This is what allows the computation 
of the ma dependence at lowest non-trivial order in perturbation theory. The actual 
order at which the computation is done however differs for valence and sea quarks. In 
the valence case, the ma dependence is computed for a free quark. For the correlation 
between two closed quark loops (fig.l.d), on the other hand, the lowest non-trivial 
order is a A s [or a 2 for a single quark loop in a background field]. This is why for 
sea quarks, we extract the ma dependence from the triangle diagram (fig. 2). So 
that in summary, given the factorization hypothesis, what we compute here for sea 
quarks is the factor replacing the (1 + ma) factor of valence quarks; the rest of 
the renormalization factor ([1 + 0(ay)] and tadpole corrections) remains unchanged. 
Another way to look at this factorization for sea quarks would be to view the gauge 
field configurations as a smooth background field with additional short wavelength 
fluctuations. In this picture, the ma dependent factor corrects the response to the 
smooth background field and the perturbative factor accounts for the short wavelength 
fluctuations. 

4 Response to continuum background fields 

The continuum computation of diagrams like the one presented in fig.2.a is rather 
straigthforward. The Feynman rules that we have used are presented in the appendix. 
These diagrams are finite and the only subtelty in their computation is the intro- 
duction of a Pauli-Villars regulator (which simply consists in substracting the same 



Kg(ma) = C(0)/ L(0,ma) 
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(9) 
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fermionic loop for an infinitely heavy quark) to ensure gauge invariance. We have 
done the computation for the axial, scalar and pseudoscalar currents (r = 7075,75 
and / respectively) which are of most direct phenomenological relevance. As men- 
tioned above, we will restrict our attention to zero momentum currents so that q = —p 
in fig. 2 . For the pseudoscalar current, however, this would lead to a vanishing result 
and we therefore give the general formula in this case. Finally, we should mention 
that when quoting our final results, we consider not only the diagram of fig. 2. a, but 
the sum of this diagram and the one with crossed photon legs [fig.2.b] (We do not 
include the external gauge fields and couplings in our expressions, since they are just 
common factors which can be brought back at the end). This gives us respectively, 
for the axial current (q = —p): 

1 p 2 

A atxv {p, -p) = ^e^apppAi— ) (10) 

7T TDj 



for the scalar current (q = —p): 



1 , s 2 \ ot P 



2 



s^ip, -p) = ^^(<W - PvP») s ^) ( n ) 

and for the pseudoscalar current: 



p»v(p, q) = ^—e»vapqaPpP{-^, \, ^-t) (12) 

&7T Tfi m z m z m z 



where the functions A,S,P are given by: 



p 



2 



if, m 2 1 , ,J—+l/2 



A ^-2\ l -V^ ln{ ^ft ]) <13) 



S(4) = 12(^/1(4) (") 

p 2 q 2 p-q f 1 [ l ~ x 2m 2 

m z m z m z Jo Jo q z x{l — x) + p z y{l — y) + 2p ■ qxy + m z 



with 




(16) 



Note that all these expressions lead to gauge invariant results (even for the axial 
current since q = —p). Finally, by using the asymptotic expansion of the functions 
A, S and P, one obtains the low momentum limits of (10), (11) and (12): 



P ->o 1 P 2 



(p, -p) — > -^2 WP/3^2" ( 17 ) 



S^{p, -p) ^ (W 2 - VpPv) (18) 

bn z m 
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P^(P,Q) P ^ ^—^apQaPp (19) 

We should mention that all the diagrams have been computed for QED rather than 
QCD. This doesn't affect the comparison between lattice and continuum at lowest 
order in perturbation theory and avoids the unnecessary burden of carrying color 
factors all the way through the computation. In fact, once the matrix element of a 
local fermionic current has been computed at lowest order in an external QED field 
and expressed in gauge invariant terms, the result can be generalized immediately to 
QCD (Color factors have to be added but there is no need to compute explicitely the 
3 or 4 external gluon lines diagram which are required to insure gauge invariance in 
QCD). 



5 Lattice results and finite ma correction factors 

We now want to compute the perturbative "response functions" on the lattice in order 
to compare them with the continuum results presented above. The computation will 
be done both for the Wilson action and for the 2-link improved action of Hamber 
and Wu (see appendix for our conventions and Feynman rules). As before, the actual 
calculations will be carried out in the context of QED. Lattice perturbation theory 
is characterized in general by the appearance of new Feynman diagrams with "multi- 
photon" vertices. In our case, the "fish diagram" (fig.2.c) appears in addition to the 
two triangle diagrams (fig. 2. a and b). The lattice expression for diagrams of the type 
of figure 2. a is: 

IIMD = J ^-Jr[TS F (k - q)T l/ (k,-q)S F (k)T II (k,p)S F (k + p)} (20) 

where the integral over k extends from — n to +n in each direction and the vertex T M 
and propagator S F are defined in the appendix. The expression for the diagram of 
fig.2.b can be obtained directly from (20) by the exchange p, /i «-> q, v. Finally, the 
fish diagram fig.2.c gives the integral 



JlAv.q) = J ^Tr{TS F (k-q)T$(k + p,-p,-q)S F (k + p)} (21) 



<Tk 

where is defined in the appendix (eq.(37)). 

Let's begin by describing the computation for the pseudoscalar current which is 
probably the simplest of all. Note that in this case (r = 75) the trace in (21) vanishes 
so that only the integral (20) has to be computed. Replacing T by 75, using the 
formulas (36) and (38) in the Appendix, and taking the Dirac trace, we get: 



d A k 1 



(2tt) 4 Z(k- q)Z(k)Z(k + p) 
{^ 5 V,{k,p)V v {k, -q)M{k - q)S p {k)S 5 {k + p) 
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-^psV^k + p)V v (k, -q)M(k)Sp(k - q)S 5 (k + p) (22) 
+^ up5 V^k + p)V u (k, -q)M(k + p)S p (k - q)S s (k) 
-4re^ a(3S V^k,p)W u (k, -q)S a (k - q)Sp(k)S s (k +p) 
-4re av/35 W^k,p)V u (k, -q)S a (k - q)Sp(k)S s (k + p)} 

The integral over d 4 k can then be carried out numerically for the Wilson or the 
improved action by using the explicit expression given in the appendix for the var- 
ious functions S,M,Z,V,W. However, since we are mostly interested in the case of 
small momenta (p,q — > 0), it is interesting to go one step further in the analytical 
evaluation and carry out the Taylor expansion of (22) up to second order in p and 
q. This computation is simplified by the fact that only the numerator needs to be 
expanded to second order (it turns out that the terms coming from the expansion 
of the denominator vanish because of the antisymmetry of the e symbol). The final 
result for Wilson fermions is (p, q — > 0): 



r d 4 k 1 / 4 

I%,(P, Q) = ^apQaPfS J ( 27r )4 Z 3 (k) (II COS fc 



v^sin 2 ^ 
M(k) -rY, 

A 



cos k\ 



(23) 



And for the 2-link improved action, we have: 

r d A k 1 ( 4 



M(k)-r^T xik) 



(24) 



where we have introduced: 



C x (k) = i(4cosA; A -cos2A; A ) (25) 
o 

T\{k) = — (4 sin k\ — 2 sin 2k\) (26) 
3 

The integrals in (23) and (24) were evaluated numerically by using a trapezoidal rule 
on a non-uniform grid. The resolution was chosen to be maximum around k = 
where the integrand peaks (especially when ma is small). In practice, we mostly used 
a set of nested (4) 4 grids where at each step, the resolution is increased by a factor of 
2 in the inner portion of the grid (see fig.3 for an example). In order to control the 
discretization errors, we further varied the resolution on each of the embedded grids 
separately and checked for convergence. The nesting procedure is stopped when the 
contribution from the inner grid becomes smaller than the overall discretization error. 
In practice, 10 levels of nesting or less were always sufficient. Once the numerical 
integrals have been computed, we can sum the contributions from the diagrams of 
fig. 2 and compare the result with the continuum formula (19). In figure 4, we plot the 
ratio of the lattice to the continuum results for the Wilson and improved action as a 
function of ma (at r = 1). All curves agree at ma = (as they should), but it is clear 
that elsewhere, Wilson fermions strongly underestimate the matrix elements of the 
pseudoscalar current. The situation is significantly improved in the case of the 2-link 
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action. These results are in agreement with those of Hamber and Wu || and Smit and 



Vink |TlJ. In tables 1 and 2, we collect the finite renormalization factors to be used in 
lattice simulations of sea quark matrix elements (inverse of the quantities plotted in 
fig. 4). We do so for the range of masses < ma < 0.2 relevant for current simulations 
and including the mass of the strange quark. The precision on these renormalization 
factors is always better than 1%. For comparison, we also include the corresponding 
factors for valence Wilson quarks (1 + ma). We would like to remind the reader that 
the numbers presented in tables 1 and 2 represent only the ma dependent part of 
the renormalization factors, the usual perturbative renormalization factors have to 
be included separetely. 

Let's now consider the case of the axial current. We have to compute (20) and (21) 
with T replaced by i"f a l5 [i is added in order to get a real answer]. Taking the trace 
directly in (20) would however lead to a rather lengthy expression. We will therefore 
restrict our attention to the case q = —p, which is sufficient for our purposes [zero 
momentum matrix elements of j a §\ . Note that for q = —p, the fish diagram for the 
axial current [(21) with T replaced by 27075] will vanish by the antisymmetry of the 
e symbol. For (20) we then obtain: 



(2tt) 4 Z 2 (k)Z(k + p) 
4{^ u5 V,(k,p)V u (k,p)[2S s (k)D(k, k+p)- S s (k + p)Z(k)] 

-2e x ^ 5 V,(k,p)V u (k,p)Sp(k + p)S s (k)S u (k) (27) 
-2e x ^ s V„(k,p)W l/ (k,p)rM(k)S^k + p)S s (k) 

where we have introduced: 

D(k, k + p) = Y J S a {k)S a {k + p) + M{k)M{k + p) (28) 

a 

As in the pseudoscalar current case, one could in principle expand the integral in 
powers of p. However this is rather difficult in the present case. One can show for 
example that the lowest non-trivial order is 0(p 3 ). Order and 2 vanish by sym- 
metrical integration and order one is a total divergence. Therefore, the expansion 
appears rather complicated and not really interesting given that the final integration 
will be computed numerically anyway. We therefore opted for the numerical compu- 
tation of (27) at various values of p and an extrapolation to p — 0. The numerical 
integration proceeds as before except that now, there are 2 points where the integrand 
peaks (k = and k = —p). We take this into account by refining the grid separately 
around these 2 points. Our integration routine generates the grids and does all the 
bookkeeping automatically. The procedure is in fact rather general and is in a sense 
a restriction of the adaptative quadrature method to the case where the position of 
"troublesome" points is known in advance. Discretization errors are controlled as 
before by doubling the resolution on each of the nested grids separately. For small 
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momentum, the results of the numerical integration can be compared with the contin- 
uum formula (17). The inverse ratio 1/ka is plotted in fig.5 and the axial correction 
factor ka presented in tables 1 and 2. The results are qualitatively similar to those 
obtained for j§ except that the lattice artifacts are smaller for the axial than for the 
pseudoscalar current. Finally, it is also interesting to look at the situation at finite pa 
as this might become useful in a more thorough analysis of lattice artifacts. That is, 
although we would hope that a numerical simulation is dominated by gluonic fields of 
small pa, we have no reason to assume this a priori and it might be important to look 
at "rough" fields too. We would therefore like to fix ma and compare the momentum 
dependence of the lattice integrals with the continuum formula (10). In fig. 6, we 
present the result of such a comparison with ma fixed at 0.2 (the continuum solid 
curve corresponds to A(p 2 /m 2 ) as given by (13) and the dots are results obtained on 
a grid of lattice momenta). We have limited the comparison to the range \pa\ < 2tt 
accessible on the lattice. 

Finally, let's consider the scalar current. Note that in this case, the "fish" diagram 
(fig.2.c) is non-zero and has to be included as well. Only the sum of all 3 diagrams in 
fig. 2 gives a gauge invariant result. In fact, we checked that the final results obtained 
in this way are proportional to a lattice transverse tensor which reduces to (5^ u — 
PuPu/p 2 ) in the low momentum limit. In our most extensive comparisons between 
lattice and continuum, however, we didn't compute the full tensor but rather focused 
our attention on a single off-diagonal element. The fish diagram, being proportional 
to 5fj, u , doesn't appear in this case and this helps in further reducing the computer 
time requirements. The expression for the triangle diagram (20) with T = 1 and 
q = —p is: 



d 4 k 



(2n)*Z 2 (k)Z(k + p) 
{V^k,p)V I/ (k,p)8M(k) [S^k + p)S u (k) + S^k)S l/ (k + p)} 
+ArV lx {k,p)W v {k,p) [2M(k)F lx (k,k+p)-S lx (k+p)Z(k)\ 
+4rW fl (k,p)V l/ (k,p)[2M(k)F u (k,k + p)-S l/ (k + p)Z(k)} (29) 
-Ar 2 W^(k, p)W„(k, p) [2M(k)E(k, k + p) + M(k + p)Z(k)} 
-45^ [V,(k,p)] 2 [2M(k)D(k, k+p)- M(k+p)Z(k)]} 

where we have introduced: 

E(k, k + p) = Y, S a {k)S a {k +p)- M{k)M{k + p) (30) 

a 

F^k, k+p) = S^(k)M(k +p) + S^(k + p)M(k) (31) 

The numerical integration of (29) proceeds as in the case of the axial current. 
The computations are made for various values of p and extrapolated down to p — 0. 
Comparing the small momentum lattice results with the continuum formula (18) gives 
the renormalization factor Ks(ma) defined in (8). The inverse of this factor is plot- 
ted in fig. 7. It appears that when ma is non-zero and moderate, Wilson fermions 
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overestimate the continuum result. The deviation is significant, being almost 30% at 
ma = 0.1. (The renormalization factor is tabulated in table 1). At large values of 
ma on the other hand, the lattice would underestimate the continuum result. The 
Hamber-Wu action provides a moderate improvement for small values of ma. How- 
ever for ma > 0.1, the situation is as bad or even worse than for Wilson fermions 
(The correction factors for the improved action are given in table 2). We also made a 
comparison of lattice and continuum results at finite momentum. By fixing ma = 0.2, 
we obtain the results presented in figure 8. The continuum solid curve corresponds 
to the function S(p 2 /m 2 ) as given by (14). Note that up to an overall factor, this is 
precisely the comparison between the functions L and C [eqs. (6) and (7)] introduced 
in section 3 [eq. (11) provides the exact correspondence]. From the curves of fig- 
ure 8, it appears that for background fields of non-zero momentum, Wilson fermions 
overestimate matrix elements of the scalar current even more than in the zero mo- 
mentum limit. It is interesting to contrast this with what happens in the case of 
naive fermions (There of course, one divides the lattice results by 16 before making 
the comparison with the continuum). The results are presented in fig. 9 and turn 
out to be much better than for the Wilson case. It is not entirely clear what the 
origin of the problem is with Wilson fermions. They probably essentially come from 
the form of the quark-gluon vertex [i.e. eq.(36)]. Gluons do not only couple to the 
quarks through 7 M but also directly through a "mass term" rW^. In fact, one can 
check that by artificially truncating the vertex (taking r = in (36)) in the Wilson 
fermion calculation of the triangle diagram and leaving the propagators unchanged, 
one obtains a result very similar to fig. 9. The result of fig. 8 indicates that one has to 
be very careful when using Wilson fermions to measure matrix elements of the scalar 
current on the lattice. Simulation results might be even more overestimated than we 
would think from the low momentum analysis given above and leading to the scalar 
correction factors given in table 1. This will have to be analyzed in more detailed 
investigations. 



6 Examples of phenomenological applications 

The results presented above indicate that the finite ma correction factors for sea 
quarks are in general rather large and will therefore strongly influence the interpreta- 
tion of simulation results. To illustrate this, we now discuss their application to two 
interesting problems [I3|. The first is the computation of the mass shift of the rj . 
In a quenched simulation, it can be estimated from a measurement of the coupling 
between the would be Goldstone bosons. The computation involves the ratio of a 
2- loop (fig.l.d) to 1-loop (fig.l.c) amplitude ||] . As a consequence of the above dis- 
cussion, we expect that they will have different renormalization factors: The 1-loop 
part behaves as a "valence" quark and will therefore pick-up a (1 + ma) correction 
for each current whereas the 2-loop diagram corresponds to a "sea" quark situation, 
in which case, the pseudoscalar factor of table 1 is relevant. In figure 10, we compare 
the mass of the rj' after correction with the raw lattice data Q . The trends are very 
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different: after correction, the mass shift of the 7/ is essentially independent of the 
mass of the quark up to the strange quark mass with a chiral limit that is in much 
better agreement with the estimate obtained from the Witten-Veneziano formula (We 
assume that the drop in the renormalized data at the lowest quark mass can be at- 
tributed to finite size effects or to the zero- mode shift effect |]ITJ). The agreement 
between the renormalized data and the Witten-Veneziano estimate Mq = 2Nfx/f% 
doesn't come as a surprise, since the correction factor that we are advocating here is 
the same as the one used by Smit and Vink [jTTJ to compute the topological suscep- 
tibility by the fermionic method [See eq.(5)]. It is just a matter of self-consistency 
that the correction factor appearing in (5) be used as well in the computation of the 
2- loop amplitude [Fig.l.d] as a function of time separation. Otherwise, the Witten- 
Veneziano relation could not be recovered on the lattice. Note also that the (quark) 
mass independence of the coupling between the would be Goldstone bosons is a basic 
assumption in the Witten-Veneziano mass formula. It is nice to be able to verify it on 
the lattice. As a second example, we mention the contribution of the strange quark to 
the mass of the nucleon. If we replace the (1 + ma) correction factor used in |5[] with 
the sea quark scalar correction factor from table 1, we find that the strange quark 
contribution to the mass of the nucleon comes down, from 30% to 21% - a change of 
30%. 



7 Conclusions 

We have computed the ma dependent part of lattice renormalization factors for sea 
quark matrix elements. The results depend strongly on the Lorentz structure of 
the current under consideration. This is unlike the valence quark situation where 
in first approximation, a single ma dependent correction factor can be used for all 
currents. We found for example that for Wilson fermions, sea quark matrix elements 
of the axial and pseudoscalar current are underestimated whereas those of the scalar 
current are overestimated. The perturbative method that we have used for computing 
the finite ma correction factors is general and could be used for other local currents, 
as well as for point-split currents. The ma dependent corrections for sea quarks 
are rather large and generally dominate the remaining part of the renormalization 
factor (2R[1 +0(av)]), which of course underlies their phenomenological importance. 
In section 6, we applied our renormalization factors to the results of recent lattice 
simulations concerning the mass shift of the rf [|J and the strange quark content 
of the nucleon ||. In both cases, it was found that the finite ma correction had a 
large impact and helped in bringing the simulation result closer to phenomenological 
estimates. 

There is also much interest in trying to reduce the magnitude of these correction 
factors by using modified actions. Our results for the axial and pseudoscalar currents 
for example give a strong motivation for using the 2-link improved action of Hamber 
and Wu in these cases. For the scalar current, on the other hand, this kind of 
improvement doesn't seem to help much and it is probably more interesting to consider 
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other alternatives, such as staggered fermions for example or working with naive 
fermions and dividing the result by 16. Finally, we would like to recall that the 
renormalization factors quoted in tables 1 and 2 have been obtained assuming that 
the lattice simulation is dominated by low momentum gauge fields. What happens 
exactly in more general situations will have to be investigated in later studies. In 
the final analysis, it seems to us that the ideas presented here would be best tested 
by doing simulations using different actions (Wilson versus improved for example or 
Wilson at different values of r), applying the appropriate correction factor in each 
case and then checking the overall compatibility of the physical results. 



8 Appendix 

The action for Wilson and two-link inproved fermions is usually defined as |9|]: 

n,fi 

(C + D ltl )Ul + ^UlM 

~ (32) 



For the Wilson action, C = D = 0; whereas for the improved action, C and D are 
chosen to cancel all terms of order {pa) 2 and (pa) 3 in the expansion of the action. 
One then finds: 

C = — D = - (33) 

It has to be remarked, however, that with this definition, the th order term and the 
terms linear in (pa) are modified by the improvement procedure. This is usually cor- 
rected at the end of the computation by introducing a wave function renormalization 
and by changing the definition of the quark mass in terms of k ( ma = ^ — 4r ). 
Here however, we will choose the alternative solution of making the correction from 
the start. That is, in the improved case, we will multiply the action (32) by 4/3. This 
will insure that our Feynman rules have the correct continuum (=low momentum) 
limit. In this way, the relation between ma and k will also be the same for the two 
actions: 

1 . . 

ma — 4r (34) 

2k 

Finally, in the same spirit of keeping things as close as possible to the continuum 
formulation, we absorb the factor of \/2k in the definition of ip. With these remarks, 



we can now introduce our Feynman rules. Our definitions are the same as in . We 
write the inverse fermion propagator in the form: 

S F 1 (p)=Y,n«S a (p) + M(p) (35) 
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the 1-photon vertex (fig. 11. a) in the form: 

F^k, q) = i~fMk, q) + rW^k, q) 
and the 2-photon vertex (fig. 11. b) in the form: 

r$(k,q,p) = V (t^VW(k,q,p) -rW<?\k,q,p 
From (35), we find for the fermion propagator: 

-iEalaS a (p) + M(p) 



Sf(p) 



Z{p) 



where we have defined 

We then have for Wilson fermions: 

S a (p) = sinp Q 
M(p) = m + 4r — r cospp 

Vfi(k, q) = cos(A; M + ^) 
W,(k,q)=sm(k, + ^) 
'1 - 

2 "' 2 



Vf\k,q,p)=sin(k, + ^ + ^ 

WW{k,q,p) =cos(^ + 
and for the improved action: 



2 2 ' 



Saip) = 7(8 sin Pa -sin2p Q ) 
6 

M(p) = m + 4r — - ^(4 cosp^ — cos 2pp) 

3 /3 



W,(k,q) = l 



4cos(/c M + y) - cos(^)cos(2A;^ + 



4 cos(/c p 



|) - 2cos(^)sin(2k, 



V™ (*, ?, p) = I [4 sin(^ + | + f)-2 cos(|) cos(^) sin(2A;, + g M + 
W< 2) (*, 9, P) = \ [4 cos(A;, + ^ + ^) - 4 cos(^) cos(^) cos(2A;, + q, + p„) 



(36) 
(37) 

(38) 
(39) 

(40) 
(41) 

(42) 
(43) 
(44) 
(45) 

(46) 
(47) 

(48) 
(49) 
(50) 
(51) 



In writing these Feynman rules, we have omitted all the factors of the lattice spacing 
a. These however are obvious on dimensional grounds and can easily be reinstated at 
the end of the computation if necessary. Finally, the continuum Feynman rules are 
obtained straightforwardly by taking the low momentum limit of the above formulas. 
The result is: 



S a (p) = p a M(p) =m V„ = 1 







(52) 
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Table Caption 

1. ma dependence of local "sea" currents renormalization factors for Wilson fer- 
mions (r=l). 

2. Same as table 1 for Hamber-Wu fermions. 



Figure Caption 

1. Quark skeleton diagrams associated with various measurements: (a) Connected 
and (b) disconnected contributions to a matrix element in the nucleon. (c) 
Connected and (d) disconnected contributions to the propagator of a flavor- 
singlet meson. 

2. Leading order diagrams contributing to the matrix element of a current in an 
external field (Diagram (c) appears only on the lattice). 

3. Example of a grid used for numerical integration. 

4. Pseudoscalar current response to a given external field. The ratio of the lattice 
and continuum responses is plotted versus ma for two choices of the lattice 
action. 

5. Same as fig.4 for the axial current. 

6. Dependence of the axial current response function on the momentum of the 
external gauge field (ma is fixed at 0.2 and pa varies between and 2ir). The 
continuum response function (solid line) is included for comparison. 

7. Same as fig.4 for the scalar current. 

8. Same as fig. 6 for the scalar current. 

9. Same as fig.8 for naive fermions (after division by 16). 

10. rj' mass shift vs quark mass after renormalization (squares) and before (plusses). 
Also included are the chiral extrapolation of the raw data (diamond) and the 
estimate from the Witten-Veneziano formula (circle), using the topological sus- 
ceptibility x obtained in ||. 

11. One and two photon vertices. 
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Table 1 



ma 


sea quarks 


pseudoscalar (ftp) 


axial (ka) 


scalar (ks) 
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